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Abstract 

This paper is concerned with the study of nonlinear stability of superposition of boundary layer 
and rarefaction wave on the two-fluid Navier-Stokes-Poisson system in the half line R+ =: (0,-|-oo). 
On account of the quasineutral assumption and the absence of the electric field for the large time 
behavior, we successfully construct the boundary layer and rarefaction wave, and then we give the 
rigorous proofs of the stability theorems on the superposition of boundary layer and rarefaction wave 
under small perturbations for the corresponding initial boundary value problem of the Navier-Stokes- 
Poisson system, only provided the strength of boundary layer is small while the strength of rarefaction 
wave can be arbitrarily large. The complexity of nonlinear composite wave leads to many complicated 
terms in the course of establishing the a priori estimates. The proofs are given by an elementary 
energy method. 
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1 Introduction 

1.1 The problem 

The dynamics of the charged particles in the collisional dusty plasma can be described by the Navier- 
Stokes-Poisson (called NSP in the sequel for simplicity) system which reads in the Eulerian coordinates 

' dtPi + dj:{p^u^) = 0 , 

Pii^^t^i “b “b dxPi^Pi^ — PiP “b PidxU-i^ 

< dtPe + dxiPeUe) = 0, ( 1 - 1 ) 

“b “b ^xPi^Pe) — PeP “b /Tedaj'Ug, 

^ dxP ~ Pi Pe- 

Here, for a = i,e, P{pa) is pressure which is given by 

P{p^)=Apl-, ( 1 . 2 ) 

where H is a positive constant and 7 ^ > 1 is the adiabatic exponent. Thus each fluid (ions or electrons) is 
regarded as an ideal polytropic gas. The unknown functions pa and Ua stand for the density and velocity 
of ions {a = i) and electrons (a = e) in plasma, respectively, and E is the electric field, while the positive 
constants p,a > 0 denote the viscosity coefficient of ions (a = i) and electrons (a = e), respectively. 
Throughout the paper, for brevity we assume 7 ^ = 7 e = 7 > 1; the case of 7 ^ ^ 7 e and 7 ^ = 7 e = 1 could 
be considered in a similar way. We also assume p-i = pe = ^ throughout the paper. One can see [I] and 
m for more information about the physical background of model (HH). 

We consider HU in the half line 1 R+ with initial data 

[p^,Ui,Pe,Ue]ix,0) = [p^o,UiO,PeO,Ueo]ix)^[p+,U+,p+,U+] aS X + 00 , (1.3) 

where > 0 and are constants. The boundary conditions are 

Ui{0,t) = Ue{0,t) = Ub < 0, Vt>0, (1.4) 

and the compatibility condition Ub = Mio(O) = Ueo(O) holds. 

In the case of Ub < 0, electrons and ions fluids flow away from the boundary {x = 0}, and thus the 
problem HU, HU) and HI in such case is called an outflow problem. The case of mj, = 0 and Ub > 0 
is called the impermeable wall problem and the inflow problem, respectively. Notice that for the inflow 
problem, there should been an additional boundary condition on the density. In the paper, we focus on 
the outflow problem in the case of Ub < 0. Here we remark that the impermeable wall problem and the 
inflow problem of the Navier-Stokes-Poisson system are left for study in the future. 

1.2 Some preliminary 

In order to study the large time behavior of solutions to the initial boundary value problem HID, HI 
and HI, we notice that in the simplified case of the electric field E = 0 and the quasineutral assumptions 
Pi = Pe and Ui = Ue for the large time behavior, the problem is reduced to consider the following single 
quasineutral Navier-Stokes equation 

{dtp + dx{,pu) = Q, 

[ p{dtu + udxu) + dxP{p) = dlu 

with initial data 

[ p , m ]( x , 0 ) = [po,uo]{x) [p+,u+], as x ^ +oo ( 1 . 6 ) 
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and the boundary condition 

u{0,t)=ub<0, Vt>0. (1.7) 

Matsumura m gave the classification of the large time behavior solutions to the outflow problem for 
Navier-Stokes equation dm in terms of (p+,u+) and Ub < 0. In what follows, let us recall some basic 
facts concerning the study of the outflow problem. The characteristic speeds of the hyperbolic part of 
(jl.5() are 


Xi=u-C{p), X 2 =u + C{p), 

where C{p) = \/P'{p) = ^/^p^ is the local sound speed. From now on, we define 


( 1 . 8 ) 


V = —, v+ = —, • • • , and so on, 
P P+ 


where v is the specific volume. Let 


c+ = C{p+) = ^p"f = 


be the sound speed and the Mach number at the far field x = +oo, respectively. The phase plane K+ x ^ 
of {v, u) can be divided into three subsets: 

^Bub ■= I (v,u) € 1R+ X R; |m| < C f- 


rtrans '■= ^ {v, u) € R+ X K; \u\ = C [ ^ 




:= < (u, u) G R+ X R; |u,| > C 


where i^sub, ^trans and flsuper are called the subsonic, transonic and supersonic regions, respectively. In 
the phase plane, we denote the curves through a right state point (ui,ui): 

BL{vi,ui) = i {v,u) € R+ X R; — = — 

1 V Vi 


i? 2 (ui,ui) = < {v,u) G R+ X R; u = ui — \Js ^ ds, v > vi 


5 '2 (ui, ui) = < (u,u) £ R+X R; u = ui + 


- P 


{yi -v), u < ui > , 


to be the boundary line, 2-rarefaction wave and 2-shock wave curves, respectively. Then the large time 
behavior of solutions to the outflow problem (IT31) . (HU) and (113 can be classified into the following four 
cases (the cases are omitted which concern shock waves): 

Case I: (u+,u+) G flsuper n{'a+ < 0} and Ub < it*. Here (Uijit*) is an intersection point of 
BL{v+,u+) and 5'2(i’+,it+), ie., 


u+ 

u+ = —u* 
v+ 



(i;+ - u*). 


_ u+ 

. 

v+ 


(1.9) 


Then there exists a unique Vb such that {vb,Ub) G BL{v+,u+), and the time asymptotic state of solution 
is a boundary layer {v,u){x) which connects (vb,Ub) with (ii+,it+), see Figure 1. By the relation of p 
and V, then we can say that boundary layer (p, u){x) connects [pb, Ub) with (p+, it+). The boundary layer 
{p,u)(x) will be explained in next section. 
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Figure 1 


Case II: (u+,tt+) G rtronsn{'^+ < 0} '^b < u+. Then there exists a unique Uf, such that 

{vb,Ub) G BL{v+,u+), and the time-asymptotic state of solution is a boundary layer {v,u){x) which 
connects (vb,Ub) with (u+,m+), see Figure 2. Here, the boundary layer (v,u)(x) is degenerate. That 
is to say boundary layer {p,u){x) connects {pb,Ub) with (p+,it+), and the boundary layer {p,u){x) is 
degenerate. 



Case III: (?;+,u+) G Hsu6p|{it+ < 0} and Ub < u+. Here is an intersection point of 

R 2 {v+,u+) and Ttrans, ie., 

t - -y + 1 ! - T"1 - T"1 

«+—vyH / s 2 ds = —\J^Av^ ^ , u^, = —\/jAv^ ^ , (1-10) 

J 

see Figure 3. This case is divided into two subcases: 

Subcase 1: If < Ub < u+, then there exists a unique Vb such that (vb^Ub) G R 2 {v+,u+), and 
the time-asymptotic state of solution is a 2-rarefaction wave which connects {vb,Ub) with 

(u+,u+), to the corresponding Riemann problem, while the 2-rarefaction wave {p ^^connects 
{pbjUb) with (p+,u+). 

Subcase 2: If m* > itb, then there exists a unique Vb such that {vb,Ub) G BL(v^,Ut), and the time- 
asymptotic state of solution is the superposition of a boundary layer (v,u)(x) connecting (vb,Ub) with 
(u*,u*), which is degenerate, and a 2-rarefaction wave connecting with {v+,u+), 

while boundary layer {p,u){x) connects {pb,Ub) with {p^,uA): and a 2-rarefaction wave 
connects {p^^^uA) with (p+,M_|_). 

Case IV: tt+ > 0 and Ub < 0. Here is an intersection point of R 2 (v+,u+) and Ttrans which 

is defined by (11.101) . see Figure 4. This case is divided into two subcases: 

Subcase 1: If u, < < 0, then there exists a unique Vb such that {vb, Ub) G R 2 {v+,u+), and the time- 

asymptotic state of solution is a 2-rarefaction wave which connects {vb, Ub) with (u+, m+). 
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to the corresponding Riemann problem, while a 2-rarefaction wave connects {pb,Ub) with 

(p+,u+). 

Subcase 2: If m* > Ub, then there exists a unique Vb such that {vb,Ub) € BL(v^,u^), and the time- 
asymptotic state of solution is the superposition of a boundary layer (v,u)(x) connecting (vb,Ub) with 
which is degenerate, and a 2-rarefaction wave connecting with {v+,u+), 

while boundary layer (p,u)(x) connects (pb,Ub) with (p^,u^), and a 2-rarefaction wave (p^^,u^^)(f) 
connects with (p+,M_|_). 



1.3 Boundary layer and rarefaction wave 

In the paper, we study the subcase 2 in Case III or Case IV without considering the other cases since 
the cases of the single wave have been studied by Duan and Yang [6] . Recalling subcase 2 in Case III or 
Case IV, there exists a unique Vb in phase plane such that (vb,Ub) G where is defined 

in (11.101). And the solution to the initial boundary value nroblem (11.11). (11.31) and dm for the outflow 
problem on two-fluid Navier-Stokes-Poisson system is expected to tend to the superposition of a degen¬ 
erate boundary layer {p,u){x) connecting {pb,Ub) with (/3*,u*) and a 2-rarefaction wave 
connecting with (p+,tt+) as t —>■ -|-c» coupling the trivial profile of electric field E = 0. 

First of all, we define the boundary layer (p, u) by the stationary solution to 

{ dx{pu)=Q, xeR+, 

pud:^u + d^P{p) = dlu, a:eM+, (l.H) 

u{Q)=Ub, (p, u)(-|-oo) = (p*,M*), inf p(a;) > 0. 


Integrating (ll.ll|) ^ over [x, -boo) for a; > 0, and letting x —?► 0, we obtain the value of p(x) at the boundary 
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{x = 0} as follows: 

Pf.= m = ^- ( 1 - 12 ) 

Ub 

Since Ub < 0, we have u^, < 0. The strength of the boundary layer {p,u){x) is measured by 

S := — Ub\. (1-13) 

In what follows let us present the existence and some known properties of the boundary layer (p, u) (x) 
connecting {pb, Ub) with (p*, it*) for the stationary problem ()1.11|1 . Here we only list the properties of the 
degenerate boundary layer. Please refer to [9] or [14] for details. 

Lemma 1.1. By the definition of (yi,,u^) in Subcase 2 in Case III or Case IV (i.e. it is located at the 
transonic curve), then there exists a solution (p, u){x) to the stationary problem (11.111) such that it = ^v, 
V = i. Moreover, u{x) is monotonically inereasing {d^u > 0) and converges to it* algebraically as x tends 
to infinity. Precisely, there exists a positive constant C such that 

(JXk+l 

\^x[p-P*^u-u^]\< fc = 0,l,2,---. (1.14) 


Since the 2-rarefaction wave [p^^,it^^] (f) is a weak solution, we shall construct a smooth approx¬ 
imation for the 2-rarefaction wave above in the following. Firstly, consider the Riemann problem for 
Burger’s equation: 

dtw wdxW = 0 , 

f W-, X < 0, (1-15) 

iii(0,a:) = wo{x) = < 

[ it;+, X > 0, 

where W- < w+. Then it is well known that ()1.15|) has a continuous weak solution (y) whose explicit 
form is given by 

W-, X < W-t, 

—, W-t < X < w+t, (1-16) 

W+, X > w+t. 



Moreover, can be approximated by the smooth function w{t,x) which is a solution to 


dtW + wdxW = 0 , 

W-, a: < 0, 

w(0,x) = wq{x) = ■^ _ 

' w- + Cq6 J y^e-ydy, x > 0, 


(1.17) 


where 5 := ic+ — ic_, g > 10 is a constant, Cq is a constant such that Cq y‘^e '^dy = 1, and e < 1 is a 
positive constant to be determined later. Then we have the following lemma. 


Lemma 1.2. Let S = w+ — W- be the wave strength of the 2-rarefaction wave. Then the problem (11.171) 
has a unique smooth solution w{x,t) which satisfies the following properties: 

(i) 0 < W- < w{x, t) < wj ^, dxW > 0 for x G M and t > 0. 

(ii) For any p (1 < p < +oo), there exists a constant Cp^q such that for t > 0 

\\dxw\\Lp < Cp,qmin{je^“p, 

< Cp^qmm{Se^~p ,6'^ 

(Hi) When x < W-t, w — W- = dxW = dfw = 0. 

(ill) lim sup |i()(x,i) — i(;^^(j)| = 0. 

i-H-oo 
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Then the smooth approximate rarefaction wave (x,t) which corresponds to the rarefaction 

wave (|) can be defined as follows: 


+ C{p'^^) = w{x, 1 + t), w- = M* + C{p^) = 0, w+ =u+ + C{p+) > 0 

1 


= u+ — /J s ^2 ds, ^ 


(1.18) 


p -2 : 


v+ = 


p+ ■ 


where w{x,t) is given in (11.171) . 

It is easy to obtain [p'^^,u^^](x,t) satisfies 

dtp^^+d,{p^^u^^)=0, 

p'^^dtvJ’^ + + d^P{p^^) = 0. 

Here we restrict [p^^ ,u'^'^]{x,t) in the half space {x > 0}. Then one has 


(1.19) 


Lemma 1.3. Let dr = \p+ — p*| + |u+ — m*| be the wave strength of the 2-rarefaction wave. Then the 
smooth approximate 2-rarefaction wave [p'^^,u^^]{x,t) constructed in (11.181) has the following properties: 
{i) dxU^^ >0,p*< p^^{x,t) < P+, M* < u^‘^{x,t) < u+, dxU^'^ ~ \dxp''^\ for x G M+ and t>0. 

(ii) For any p {1 < p < +oo), there exists a constant Cp^q such that for t > 0, 


|| 5 x[p’'^u’'^]||LP(R+) < C'p,gmin{5re^ p,<5r(l + t) 




Wdlip"'^, u''^]\\lp(r+) < Cp^qmm{Sr€^ p,6.. 


2-J 


F-- + 7 


(l+t)-^+i}. 


(in) [p’’2,-«’’2](0,t) = [p*,it*]. 

(iv) lim sup |[p''^,'u’’^](x,t) — \p^^,u^^] (7)! =0. 
t^+^xeR+ 

Now, we define 

[P,u]{x,t) := [p,'u](x) + [p’'",'u’'^](x,t) - [p*,u,]. 

By a straightforward calculation, we have 

dtp + dx{pu) = /, (x, t) e R+ X R+, 
p{dtu + udxii) + dxP{p) = dlu + p, (x, t) e M+ x R+, 
(p,u)(x,0)-)■ (p+,u+), as X-)■+00, (p, u)(0,t) = (p6,Mb). 


( 1 . 20 ) 


( 1 . 21 ) 


where 


/ =dxp{u^^ - M*) + dxu{p"‘" - p*) + dxp"'^ {u - u^) + dxu"^^ (p - p*), 
g = — + udxu{p'~^ — p*) + p [dxu{u^^ — u,) + dxu’"^{u — u*)] 

+ dxP [P'{p) - P'(p)] + dxP^^ [P'{p) - P'{p^^)] - ^^dxP^HP - P*)- 


( 1 . 22 ) 


From (11.111) ^ and (11.201) . it is easy to know 


l/l + Ip + d1u'"'^\ < C{dxuipP^ - u*) + dxU^^{u^: - u)} , 

|a./| < C {{\dlu\ + (dxhfXu^^ - u,) + dxudxu^^ + \dyX + {dxU^X^} 


where dxU > 0, dxU^^ > 0 and u < it* < it’’^. 


(1.23) 











1.4 Main results 


We can easily derive E{x,t) = — J^^[pi{y,t) — Pe{y,t)]dy from (II.IL if we assume that E{x,t) —>• 
0 as a; —>■ +oo holds. Then we can define E{x, 0) = — J^°°[pio{,y) — Peo{y)]dy- Now we are in a position 
to state our main results. 

Theorem 1.1. Let a = i,e and assume that constant states Ub, u* and the infinite state (p+,u+) satisfy 
Subcase 2 either in Case III or in Case IV. There exist some positive constants eo > 0 and Cq > 0 such 
that if 


||[p„o(-) - p(-,0),u.o(-) - «(-,0)]||^, + !!£;(., 0)f + +55 < el (1.24) 


where e > 0 is the parameter appearing in (I1.17|) . then the initial boundary value problem (imi . (foi) and 
dn admits a unique global solution [pa,Ua, E]{xC) satisfying 


S^p\\[pa - P,Ua - U,E]{-,t)\\fji <CoSo. ( 125 ) 

t>0 ' ■ ' 


Moreover, the 


solution [pa,u 


lim sup 

t->- + 00 


a,E](x,t) tends time-asymptotically to the composite 
[pa,Ua\{x,t) - [p,u] (x) - [p^Cu^l 


wave in the sense that 
= 0, (1.26) 


and 


lim sup \E\ = 0 . /I 97'i 

As it is well known that, there have been a great number of mathematical studies about the outflow 
problem, impermeable wall problem and inflow problem of Navier-Stokes system, please referring to 
[3 [U H [ini [13 [E] and the references therein. However, to the best of our knowledge, there are very 
few results about the above mentioned problems for NSP system. Duan-Yang [3 firstly proved the 
stability of rarefaction wave and boundary layer for outflow problem on the two-fluid NSP system. One 
important point used in (3 is that the large time behavior of the electric fields is trivial and hence 
the two fluids indeed have the same asymptotic profiles which are constructed from the Navier-Stokes 
equations without any force under the assumptions that all physical parameters in the model must be unit, 
which is obviously impractical since ions and electrons generally have different masses. The convergence 
rate of corresponding solutions toward the stationary solution was obtained by Zhou-Li [20) . In the 
paper, we study the nonlinear stability of the superposition of boundary layer and rarefaction wave for 
outflow problem on two-fluid NSP system. The complexity of nonlinear composite wave leads to many 
complicated terms in the course of establishing the a priori estimates. Lemma 14.11 plays crucial role to 
deal with the complicated terms. Compared with Navier-Stokes system, the key to prove Theorem ll.il 
for NSP system is to deal with the extra electric field E which is no longer integrate in space and 
time due to the structure of the Poisson equation in (12.2L . The detailed way to deal with the terms 
involved with electric field E is stated in (I2.12I) . (I2.13D and (12.151) . Finally, we remark that NSP system 
CH) in the non-dimensional form depends generally on the ratios of masses, charges and temperatures 
of two fluids. If we don’t ignore these physical coefficients, the two-fluid plasma system exhibits more 
complex coupling structure and the corresponding analysis of the large time behavior of solutions becomes 
more complicated, referring to [3] and (3- Hence it is meaningful and interesting to study the general 
physical situation for the nonlinear stability of superposition of boundary layer and rarefaction wave on 
the two-fluid NSP system in the future. 

Finally, we refer readers to [3I3I3II3I13 and references therein for the study of the related works on 
the NSP system. Here we would still mention several most closely related papers: [nmi] for the spectral 
analysis and time-decay of the NSP system around the constant states, mm for the global existence 
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of strong solutions to the one-dimensional NSP system with large data. Recently, the stability of the 
superposition of rarefaction wave and contact discontinuity for the NSP system with free boundary has 
been obtained by Ruan-Yin-Zhu m- For the investigations in the stability of the rarefaction wave of 
the related models, see also [4] for the study of the more complicated Vlasov-Poisson-Boltzmann system. 

The rest of the paper is arranged as follows. In the main part Section 2, we give the a priori estimates 
on the solutions of the perturbative equations. The structure of Poisson equation and the symmetry of 
two-fluid system play important roles in the proof of the a priori estimates. The proof of Theorem ll.ll is 
concluded in Section 3. 

Notations: Throughout this paper, C denotes some positive constant (generally large) and c denotes 
some positive constant (generally small), where both C and c may take different values in different places. 
jj) _ 2 jP(]R_|_) (1 < p < -boo) denotes the usual Lebesgue space on R+ with its norm || • ||lp, and when 
p = 2, -boo, we write || • ||l 2 (r^) = || • || and || • ||l<-(r+) = 1| • ||oo- We use = Ff®(K+) (s > 0) to denote 
the usual Sobolev space with respect to x variable. 


2 The proof of a priori estimates 

Let [pi^UiT pe,Ue,E] be the solution of the one-dimensional two-fluid Navier-Stokes-Poisson system (II.ID . 
dni) and dm). Let [p, m] be the solution of p.2ip . Now, we put the perturbation ■i/ij, pe, ■^e] by 


fi= Pi- P, il^=Ui- U, ip^= p^- p, -06 = Me - u- 

Then, from dm and (|1.21|) . [pi, 0,, pe, V'e] satisfies 

dtPi + Uidxifii + Pid^'ipi = -fi, 

Pii^t'^i “b Uidx'^pi^ P {pi^dx^i — 9i P PiE i 

dtPe + UedxiPe + Pedxfpe = “/e, 

Pe{,^t'4^e “b Me(93,0e) -b P (Pe)^xb^e ^x'4^e 9e PeE^ 
dxE = ipi- ipe, X e R+, t > 0, 

(0i,0e)(O,t) = 0, 

_ (0i,0i,0e, V'e)(a;,0)-)■ 0, as a;-)■-boo, 
where fa, Pa [ot = i, e) are the nonlinear terms, given by 

fa = dxilVa + dxPIpa + /, 

9a = PadxUlpa + OxP [P'(Pa) “ ^’'(p)] + “ dxPip)]‘^ + 9 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


We define the solution space X{0,T) by 

Y(0,r) :={[ipa,f^a,E] € C([0,T];iLi), [dx^a,dxE] € 

dxifa & 0«(o,t) = o, a = i,e, V(a:, t) £ [0,-boo} x [0, T]} . 

The local existence of (12.2D can be established by the standard iteration argument and hence will be 
skipped in the paper. To obtain the global existence part of Theorem ll.il it suffices to prove the following 
Proposition 12.II (a priori estimates). 

Proposition 2.1. (a priori estimates). Assume all the conditions listed in Theorem 11.11 hold. Let 
[ipi,'>pi,ipe,fje, E] be a solution to the initial boundary value problem (j2.2|) onO <t <T for some positive 
constant T. There exist some positive constants C and £i such that if 


sup (||[(pi,0i,(pe,0< 

0<i<T 


^1 + llif(t)ll) -b e -b d < El, 


(2.4) 







10 


then the solution ipe,ipe, E] satisfies 


sup \\[ipi,'ipi,(pe,'fe,E]\\]ji + \\^/dxU[lfi,'ipi,(Pe,'(pe]\fdt+ / \\da;[(pi, ipe, E]\\'^ + \\dx [tfi, tpe] Wnidt 


0<t<T 


^11 [<fi0,tpi0^ PeO, V'eo]||ffi + ll-E’(O) Oil ^ 59 ^ . 

Using (12.41) and the following Sobolev inequality 

|h(a;)| < V2\\h\\^\\hj:\\^ for h{x) € H'^{R+), 

we have 

\\[(pi,1pi, ipe,-felWoo < \/2ei, 


(2.5) 


( 2 . 6 ) 


(2.7) 


which will be repeatedly used in the following. 

We prove Proposition l2.1l bv elementary energy methods. Lemma |4.II in the appendix plays a key role 
in the stability analysis. Before our estimates, we should point out that the general constant C below 
may depend on the strength of the rarefaction wave 6r since the rarefaction wave considered here is not 
weak. Now, we prove Proposition 12.1 1 bv the following three steps. 

Stepl: The zero-order energy estimates. 

For a = i, e, we define the function 


= 


^’(PcP) = [ 
J p 


P[s) - P{p) 


ds 


and rja = Pa^a + ^PafJa- Dhect calculations give rise to 

dtVi + dx [uiPi + {P{pf) - P{p))'4’i - ■fxdx'fi] + dxU [P{pi) - P{p) - P'{p)>Pz + Pi'f'^] 

I /Q / \2 I rp ^ Pif’i u'! 

+ [Ox'lpi) =Pi 1 p^E-^^U— - g— - P{p)f — 

P P P 


( 2 . 8 ) 


and 


dtVe + dx [UePe + {P{Pe) “ P{p))lpe “ i’edx'ife] + OxU [P(pe) “ P(/3) “ P'{p)Pe + Pet/’e] 

I /n / \2 ; 771 a2 - ‘PeV'e ~ PeV’e 

+ KOxVe) = -PelpeE - d^U— - g— - P {p) f — . 

P P P 


(2.9) 


Taking the summation of (12.81) and (12.91) . and integrating the resulting equation with respect to x 
over R+, we arrive at 

4 : / (rii + ge)dx + \ub\[{pi^i)lS>,f) + {pe^f){Q,t)]+ [ [{dx'fi)'^ + {Oxipe^] dx 

dt Js_+ JR+ 

+ f dxU [P{pi) - P{p) - P'{p)pz + P{pe) - Pip) - P'ip)Pe + Pi'ff + PeV'e] dx 
J ]R_|_ 

= [ [Pi'^i - Pe4’e)Edx - f (dx 

Js.+ Jr+ \ p p J 


h 

.Pi'f’i , ^Pe'ft 

g^-+9—r- 


dx- I [P'{p)f)^ + P\p)fi£l]dx. 


( 2 . 10 ) 


Q 2 


Q 3 


Here we have used the boundary condition (12.2|) ^ and Uf, < 0. 
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From Poisson equation and mass conservation equation, we have 


d:cE = Pi-pe, dtE = PeUe - PiUi. 


( 2 . 11 ) 


Now we mainly make use of (12.111) to deal with the difficult term Ii. Then one has by integration by 
parts 


= 


/ E{piUi - PeUejdx - / E{pi - Pe)udx 

JR, Jr. 


1 d 




dr.uE^dx. 


I 2 


( 2 . 12 ) 


Notice that dxU = dxU + dxU^'^ > 0 from dxU > 0 and dxE''^ > 0. Now we pay our attention on the 
bad term I 2 since the electric field E is no longer integrate in space and time due to the structure of 
the Poisson equation. The main idea is to make use of the good term 



dxU [pi^^i + Pei’l] dx 


to absorb l 2 - For this, multiplying (12.2L and (12.21) ^ by ^EdxU and —^EdxU respectively, then 
integrating the resulting equations over ]R_|_ and taking the summation of the resulting equations, one has 


I 2 = 


l_d 
4 dt 


[ dxu{^i-il^e)Edx-^ [ {'i(ji-'il;e)dtEdxudx-^ ( {'ipi - tpe) Edtdxudx 
7 m , 4 4 


f '^e^x'4^e') EOxUdx f 

4 Jr, 4 


^3 

1 f (dl'ipi d^tp. 


Pi Pe 


Edxudx 


h 


+ - 


1 f P'iPe) 


dx^i - -dx^e Edxudx + 


1 


Pi Pe 


Edxudx. 


(2.13) 


^8 


From (I2.11L we have 


dtE = Petpe - Pii>i + {^e “ ^i)u. (2-14) 

Then we make use of (12.141) to deal with the difficult term I 3 . Therefore, one has 

-^3=7 / dxuipi'll^l ^ Pe'^l)dx-\ I dxu{pe Pi)'lpi'ipedx[ dxUi'lpe - - ^i)udx . (2.15) 

4 7 m^ 4 4 


h ho 

Combining (I2.10I) - (I2.15|) . we arrive at the following equality 


_d 

dt 


Id 

- Ux-J- 


Ill 


E^ 

Pi$i(0,t) + pe$e(0,t) + ^(0,t) 


dxuii’i - 4’e)Edx + \ub 

+ [ [{dxipif + {dx'ipef]dx + f dxu[Pipi) - p{p) - P'ip)y:>i + Pipe) - Pip) - P'ip)^e]dx 

J m_l J m_i_ 


/ dxU [pilpf + Pelpl] dx - Ig- ho 

Jr^ 


— Qi d- Ii + /ll. 


(2.16) 


2=1 2=4 
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First of all, we use dSH) to deal with the left terms in (I2.16|l as follows: 


\Ub\ 


+ ped>e(0,t) + ^(0,t) 




( dj,U [P[pi) - P{p) - P'{p)<Pi + P{Pe) - P{p) - P'{p)(fe] dx > c\\^/^^[pi,^pe]\\'^, 


and 


3 1 3 ' 

+^{Pe+ Pi) V'iV'e + ■^Pei’l 


> 


> 


/ da:U [pill)} + Pelpl] dx - Ig - ho = d, 

Jr+ J'r+ 

j [ pd^u[3lpf +2'lljitl}e+3hl]dx- C\\[pi,Pe]\\oc,\\\/^['(pi,he]\\ 
4 7k+ 

i f pdxu[2{tlJi +hl) + {hi+'Pe)'^]dx-C£i\\\/^[tpi,'lpeW- 


dx 


Therefore, we have 


d^^il [pill^i + Pei’l] dx- Ig- ho > 


where we take ei small enough. 

Before our estimates, we take g = 10 and 0 = i in the following for brevity. By employing (12.71) . (12.41) . 
(I 2 . 3 [) 2 . (12.2L . Lemma lOl Lemma HTTl Young inequality, Cauchy-Schwarz’s inequality with 0 < 77 < 1, 
Sobolev inequality (12.61) . the boundary condition = V'e(0,t) = 0 and integrating by parts, we 

obtain the estimates on the right terms in (12.161) as follows: 

IQil + IQ2I + IQ3I 

<C\\[pi,Pe,^i,^e]\\oo [ (\f\ + \g + dlu^^\ + \dlu^^\)dx + C [ \dlu\ + Ifl + hi + he) dx 


<(711 [pi, Pe, hi, he] II = \\dx[pi,Pe, hi, '*/'e]|l ^ 


1 + St 


+ e®(l + ln(l + k) + (1 + 


+ ck [p^io,t) + pI{o, t)] + C5\\dh^,, pe,hi,he]f 


-10 
0 9 


<C{k + e^°)\\dx[pi, Pe, hi, 4’e]k + C'—-+ Ce^o (1 pt) + Ck \p1{0, t) + pliO, t)] 

(1 + dt )3 


\h\ <C||[ 7 / 7 ,,V 7 e]||oo||i?||oo||[^^t 9 .M'-=]||L 7 

<CeHl + t)-^+-h\[hi,he]\\-hmh^,M\-^\\E\\i\^^^^^ 

<Ce^{l + t)-i +Ce^\\^x[h^,he,E]f, 


llo] + l^el + \h 

<C\\dxU'"'^\\oo\\E\\\\dx[hi,he, dxhi, dxhe, Vi, V’eWl 

+ E \dx[hi,h ei hi,dxhe,Pi,Pe]\ \E\dxUdx 
«/R_|. 

<(7e®(l + t)~^'^~^'>\\E\\\\dx[hi,he,dxhi,dxhe,‘Pi,‘Pe]\\ 

+ ckdx[hi, he, dxhi, dxhe, 'Pi, Pe, E] |p + CkE'^{Q, t) 

<C{5 + e^)\\dx[hi,he,dxhi,dxhe,Pi,Pe,E]k + CkE‘^{{)h) + Ce^{l+t)~i, 
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|4| <C / {\gi\ + \ge\) \E\\dxu\dx 


and 


<C / {\da;U\{\tpi\ + \lpe\ + I'fil + \<fie\) + \dlu\{\(pi\ + \(pe\)} \E\\da:U\dx + C / \g\\E\\da;U\dx 

<c / {\d,u\^ + \dlu\ + + \dlvr^\) {iPl + iPl +^^ + ^l+E^)dx + C [ \g\^dx 

Jr+ Jr+ 

<CP [iflio, t) + (^2(0, t) + E‘^{0, t)] + C5\\d^[<^i, i’e, E]f + Ce^+t (1 + + C5{1 + 1 )-^ 

+ C'(||[v5i,(/3e,V'i,'*/'e,£']|| 119:1, [(/5i,(/3e,V'i,'*/'e,£']||) (||9^w’'" ||il + ||9a;'u’'"f) 

<C {6 + eio )||93;[(^i, ipe, ipi, ipe, E]\\^ + C(eio + (5)(1 + t) ^ -(- C5^ t) + p^{0, t) + E^(0, i)] 

1 -^ 11 1 / \d:^u\{\lpi\ + \l{je\)\da;E\dx + C [ \d^u''^\ {Itpil + \lpe\) \da:E\dx 

<g\\da;E\\'^ +Cr, [ |9xMp + IV'en + C'||9xM’'^||oo(||'*/'i|| + II II ) II 9x^^|| 

Jr+ 

<g\\d,Ef + ^PeW + Ce^il + <)-'^'-"^(IIV'*f + M) + Ce^\\d,Ef 

<(?7 + C'e8)||9xi?|p + CrjS^\\dx['4’i, '0e]||^ + C’e8(l+t) 

Substituting the estimates above into (I2.16|) and integrating the resulting inequality over [0,T] and 
using Cauchy Schwarz’s inequality, and taking e, S and Si small enough, one can see that 


rT , 


dt 


(2.17) 


\\[pi,Pe,'ipi,-ipe,E]f + / \\dx['tpi,1peW + \\\/^[Pi,Pe,'ipi,-ipe]\\ 

Jo *- 

+ [ + i‘Pe)'^{0,t) + E'^{0,t)]dt 

Jo 

<C (||[(pio,‘/5eo,l/’*o,'0eo]||^ + ||^^(a;,0)||^) + (r? + Cera + C 6 i)\\dx[pi,Pe,dxipi,dxipe,E]\\'^ + C(era + S^). 
Step 2. Dissipation of dx[~Pi, ‘Pe:E]. 

We first differentiate (12.21) ^ and (ESDa with respect to x, respectively, to obtain 
dtdxPi + dxUidxPi+Uidlpi + dxPidx'tpi +Pidl'ipi + dl'upi+dxudxPi + dxpdx'ipi+ 81^4)1 + 8 x 1 = 0 (2.18) 
and 

8t8xPe + 8x'ae8xPe+^e8j.Pe+8xPe8x4e+ Pe8x4e + 8xdl'Pe + 8xn8xPe + 8xp8x4e+8xP4e. + 8xf — 0- (2-19) 
Then multiplying (12.21) ^. (12.21) .,. (12.21) .. (12.181) and (12.191) by 8xE, d^vB ^ inte- 

O ^ ^ pi pg pg 

grating the resulting equalities over R+, one has 

f {8xEYdx = [(^e(0, t) - piif), t)]E{0, t)- ( 8x{pt - Pe)Edx, 

Jr+ Jr+ 

^ ^2 


8t4i8xPidx + f Ui8x4i8xPidx + f ^ {8xPi)^dx 
Jr^ Jr+ Pi 

= [ dxPiEdx + [ 8l4i—^dx— [ gi^di^dx. 




8t4e8xPedx + I Ue8x4e8xPedx + 


/r+ P 
P'iPe) 


Pe 


i8xPe) dx 


8xPeEdx + ( 8x4>e—^dx— f ge-d^dx, 

Jr+ Pe Jr, Pe 
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^dtd,^,dx+ [ d,uS^^dx+ [ u?-^^^^^S^dx 

K+ PI Jb+ pi Jr+ Pi 

+ [ ^^^dxPidxi’idx + f dx 

JR+ Pi Jr+ Pi 

f ^2 / j f 02'' 9x^i j f a -o / 

= - ^x'^i -- —dx- / OxPOxWi -^ 

•/R-I- Pi *1 r+ Pi *lR-f- Pi 

^2 ^x’pi f ^ x^x^i 


dx 


and 


dipi>^^dx- dxf^dx 

Pi Jr+ Pi 

9x‘fier,^ , f . (5x<Pe)^ 


f ^^^dtdx(Pedx + [ dxUe dx + f U, 

Jm.^ Pe ^M-u Pe */R-i- 


K+ Pe •/R+ P( 

dxPe a _ c\ I j_„ if c\ '• i^xPe^ 


dx‘fedl‘fe 

pI 


dx 


+ 


( n 2 

^R+ Pe 


f ( 

~dxPedx'4^edx ~h / 

d R-i- 


da; 


+ P~ 

d^ipe-^^dx - I d^uife-^P^dx - I dxpdxip. 


,2 / P^x^e 1 f 02 ''. ,_ j™ f n „/. ^x^e 


dx 


f dlpipe^^^dx - j dxf^^^dx. 

R_i. Pe «/lR_i_ Pe 


x_|. re «/M+ 

The summation of the equalities above further implies 


j-i [ (pPidxPi “t“ '^edxPe') dx T , ^ ( n 2 n 2 (^a:;‘Pe) ) dx 

dtJui dt \ 2 p^ 2 pi ) 


dx'd 


{dx<Pif , , P'{Pi),f, x2 , P'{Pe),r, '2 , /o z7\2 


+ 


Pe Pt Pe 


dx 


= [V3e(0,i) - I^i(0,i)]£’(0,l) + [ {ipidtdx(Pi + ipedtdx(Pe)dx 

Js.+ 

/ Pi P^tPi “t“ (^x^e) Pe ^tPe) dx / {Uidx‘4^i^x^i ~\~ da; 

aRi Jr, 


dx^i dx^e ^ , 

Pi-h 5e- I dx - 

+ \ Pi Pe J ,/ M+ 

dx^idl'Pi , dx^edl^. 


Pt 


+ Me 


dxUi 


dx — 


i^dx^pi) 


~\~ dxUi 


{dxPeY 


dx 


Pi Pe 

9xPi o Q ; I ^xPe o a ; 1 J 

2 ^xPi^xPi 4” 2 ^xPeOxPe J dx 


R+ V Pi 


,2-' P^xPi o2'',,„ ^xPe\ f f xi -o „/. P^xPi , o ^o „/. ^xPe 


d 2 ^ '^x'ri I o 2 ^ 

p* pf 


dx- [dxpdxtIJi - — + dxpdxil)e 2 

P* Pe 


dx 


o2 - ; '9a;<Pi o2 - ; (^xPe \ , 

OxPVi^r + ^xPVe -5- dx - 

+ V Pi Pe 


12 


+ dx = J 2 Ji, 

Pi Pe J 


( 2 . 20 ) 


where Ji {1 < I < 12 ) denote the corresponding terms on the left hand side of ( 12 . 201 ) . 

Notice the fact that \dxn\ < CdxU, Id^n] < C(|i9^m| + \dxu\‘^) and ut < 0. We now turn to estimate J; 
(1 < Z < 12) term by term. By applying Holder inequality, Cauchy-Schwarz’s inequality with 0 < 77 < 1, 
Sobolev inequality (12.61) . Lemma [T31 Lemma ITTl ( 12 .4L (12.71) . (11.11) ^. (II.ll) ^- ( 11 . 211 ) ^. (I1.14L (12.31) . the 
boundary condition ipi{Q,t) = ipe{^,t) = 0 , and integrating by parts, it is direct to derive the following 
estimates: 


|di|<^?( 0 ,t) + ;^^( 0 ,t) + L; 2 ( 0 ,t), 
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J 2 = dxi^idxipiUi - pu)dx / dxi)edx{Pe'^e- Pu)dx - / +'0e)9a;/(ix 

JE 4 . JE 4 . 7]R+ 

= / + / dxpii^idxi^i + 'llJedx'lpe)dx 

Jr+ Jr+ 

I [(^xV^z) ^2 H” (^xV^e) “1“ / {,Pi^x^^x'4^i “t“ PeQxU0x'4^e}dx 

Jr+ Jr+ 

“t“ / ip^idx'4^i^xPi ~\~ '^e^x'4^e^xPe}dx / H“ '4^e')9xf dx 

Jr+ Jr+ 

<V\\dx[‘Pi^‘PeW + Cr,\\d^['lpi,'tl}eW + C{e + 6)\\^/^[(pi,(pe,1pi,1pe]f + C\\[tjji,1pe]\\oa\\dxf\\L^ 

<V\\dxyi,^eW + Cn\\dx['>pi,'ll^eW + + ^)|| (/£>e, V'i, V'e] IP 

+ C||[l/^i,'!/'e]||^||9x[lAi,'0e]||^ ^(1 + ty^ + e®(l + + e?(l + 

<»7ll^x[V5i,</?e]|P + Cr,\\d^[tpiyeW + C{e + 5) || , (/Jg, , ^e] |1 ^ + 0(6^ + eA)(l + t)~i, 


>^3 + >^6 + J 7 
1 




< 


P?( 0 ,i) 

1 


p 2 ( 0 ,<) 


^dxiiid^y^p, ^ + -dxu{d:c^eYPe ] dx + 


“1“ / ( ^px^^ii^x^i) Pi “t" 2^3;^e(^x^e) Pe ) 


^^x'4^iidx^i) Pi ~t“ pdx'^ei^x^e) Pe ) 


<C'(e + ^)l|5x[</?i,</5e]|P + C'||a;:c[V'i,V'e]ll"|19^[V'i,V'e]ll"l|9xbi,</5e]f 

<(7(6 + J + £l)(||aa:[(/7i, V?e]f ) + C'ei||a^[V’i, f/^elf, 

I J 4 I + |7lo| < (?? + C5 + C'e)||9x[‘/^i; ¥^e]|p + {Crj + CS + Ce)\\dx['tpi, '!/^e]|Pj 

\J5\ <c [ {\dlu\ + \dxp\ + dxu)\y,^ e; V’i)'0e]||9x[v?i, i^e]Ma; + C\\dxy,'Pe 


<p\\dxVPi,‘fe]\\ +<7^(11311 +\\dxUy^ + \\dxvJ'^\\yj+C5\\dx[^Pu‘Peyiye]\\ +C5 y+^ey,t) 
<(77 + C5)||9a;[(/3i, (/3e,'i/’i, ■*/'e]|P + i) + (/?g(0, <)] + {1 + t) ‘'^+5(l + <) ^ 

yv + C6)\\dx[<Pi,'Pe,AyeW + CS'^[(pf{0,t) + y{0,t)] + cy + ei){l+t)~i, 

l^sl <C\\dxp\\oo\\dxy,‘Pe]\\\\dx[tpiye]\\ + [V^i, V'e] II ^ [V'i, V'e] II ^ , ¥>e] f 

<(7(e + (5 + ei)\\dxy,‘Pe, V'i, V’elf + C'SlIIS^IV'i, V'elf, 

I^qI + I Jill <(7 / (|a^-«| + |9a:Mp)|[(^i,(^e,V'i,V'e]l|5xbz,‘/5e]|fia; 

+ C f [\dxU'^^\ + \dxU’^^P\[(pi,ipeyi,'lpe]\\dx[‘Pi,<Pe]\dx 
Jr+ 

<cs'^[y{o,t) + <fi{o,t)] + c^dx[(Pi,(peyiyeW 


+ c 


:^'''‘*(1+^) + e^^(l + t) II I'i/^i,'0e]|| ||9x[</?z, V’e 


and 


<(7(^ + eO||93:[(/3i,(/5e,-04, V’elll^ +<76^(1 +t) =, 


IJ 12 I <||9x[(/5z,¥5e]||||5^/|| < p\\dx[(pi,>PeW + CrjWdxfW" 

^- 2 (l-i) 


—''711^2: [^ii ¥^e] II + C'l^e ’(1+^) ® + Cr;((5 + e)(l + i) 

^VWdxi^PiT^Pey + Cj^{6 + €){1 + t) 
where we take g = 10 and 0 = 7 in the above estimates. 


-2 
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Inserting the above estimations for Ji {1 < I < 12) into (12.2011 and then choosing ei, e, S and 77 so 
small, and integrating (j2.20l) over [0,T] and using (12.171) . Cauchy-Schwarz’s inequality with 0 < 77 < 1, 
one can see that 


\\dj;[ipi,(fe]f + [ \\\/^dj;[ipi,(pe]\\'^dt+ [ \\d:^[ipi,ipe,E]\\'^dt 


^||['0zOj V'eo]|| + \\E{x, 0)|p + II [lyJiOj '4^eo \||^i j + {rj + + C6^ + ei)||9^[7/>i, T/je] Ip + (7(6 1“ + ^ 9 ). 

( 2 . 21 ) 

Step 3. Dissipation of d'^[ipi,'ipe]- 

Multiplying (l2.2D o and (12.21) ^ by —^ and —respectively, and then integrating the resulting 
equations over R_|_ and taking the summation of the resulting equations, one has 


d 

dt 


+ 


[ dx+ f 

\ / t/M.. 


\ p^ 


Pe 


dx 


Edl{ipi — il)e)dx + 


P'{pi) a a2 ; I P'{Pe) ^ a2 ; 1 j 


P^ 


Pe 


[ {uAipidlip^ + Uedooipedlipe) dx + f (—dlf)i + — dlf)^ dx 

JK.+ \Pi Pe J 

16 


( 2 . 22 ) 


i=13 


where we have used the boundary condition '0^(0, t) = 7/’e(0,t) = 0. 

We now turn to estimate J; (13 < J < 16) term by term. By applying Cauchy-Schwarz’s inequality 
with 0 < 77 < 1, Sobolev inequality (12.6L Lemma [1.31 Lemma [4.11 (11.1411 and integrating by parts, we 
can obtain that 


Ji3 =E{Q,t)[{d^ij)i){{),t) - {dxij)^){Q,t)] + f d^Ed^ipfi - V'e)dx 

<p[{dx'>fif{Q,t) + + C,,L;^(0,t) -b ]^\\dx[f)i,1peW + \\\9xE\\‘^ 

<??l|9x[l/’i,V’e]llL + C^E^{0,t) + ]^\\dx[i)i,ifeW + \\\dxE\\'^ 

<ri{\\dx[^M\^ + \\dl[f;M?) + Cx,E\t),t) + + \\\dxE\\\ 

\Jl4\ + IJlsI < v\\dl[f;i,'ll)eW + Cr,\\dx[^Pi,^e,'^i,f^eW 


l^iel < ( 7 ? -b C5)\\dx[(pi,‘Pe,'>l’i,'4^e,dxii)i,dx'>fe]\\'^ + C6'^[(pi{Q,t) -b ipl{0,t)] + Cjj{S + e^{l + t) ^ 

where we take g = 10 in the above estimates and the estimate of Jie is the same as J5. 

Inserting the above estimations for J; (13 < J < 16) into (I2.22p and then integrating (12.221) over [0, T] 
and using (12.171) and (I2.21L one can see that 

\\dx[f}i,llJe]f + [ \\dl[llji,fjeWdt 

Jo 

^11 [(^iO, <7260,'0iO,'0eo] lllfi + II 0)1 II ^ -b C ^£10 -b ^9^ . (2.23) 

where we choose ei, e, S and 77 sufficiently small. 
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Proof of Proposition \Kl[ Now, following Step 1, Step 2 and Step 3, we are ready to prove Proposition 
12.1[ Summing up the estimates (|2.17p . (|2.2ip . (|2.23p and taking e, S, r] suitably small, we have 

pT 

sup (\\[ipi,ipe,tlJi,1pe]\\m + \\Ef)+ \\y/K^[ipi,(pe,'ipi,'4’eWdt 

0<t<T ^ ^ Jo 

+ [ \\d4<fi,'fe,E]\\‘^dt+ f \\dx[tpi,tpe]\\Hidt 

Jo Jo 

— ^ (l|[V5iO, ^e 0 , 1 pi 0 , V’eolll^i + || [E{x, 0)]||^ + £1“ + ^9^ . (2.24) 

From (12.2L . it follows 

\\d.Ef<\\[^.,^e]f, 

this and (I2.24|) imply the desired estimate (12.51) . Thus the proof of Proposition [23 is completed. □ 


3 Global existence and large time behavior 

We are now in a position to complete the proof of Theorem 11.11 

Proof of Theorem 11.11 By the a priori estimates (1^ . there exists a positive constant Cq such that 

\\[ipi,iPe,if^,ipe,E]fj^i < Co (\\[ipio,Teo,'ipio,'ipeo]\\Hi + ||[£i(a;,0)]|1 ^ + £w (3.1) 

holds. It is straightforward to see that there exists a small constant Eq such that if 

\\[ip^o,Teo,tp^o,tpeo]\\Hl + \\E{0,x)f < eg, 


we can close the a priori assumption (12.41) by choosing ei = 4^Co(eg + £io + ds). By letting e and S be 


small enough, then the global existence of the solution of p.2|) follows from the standard continuation 
argument based on the local existence and the a priori estimates in Proposition 12.II Moreover, dSlD and 
(11.241) imply (11.251) . Our intention next is to prove the large time behavior as (11.261) and (11.271) . For this, 
we first justify the following limits: 

|2 


and 


lim \\dj:[(pi,(pe,'lpi,i^e]it )\\^2 =0, 

t^ + OO 


lim \\d^E{t)\\ =0. 

t^ + OO 


(3.2) 

(3.3) 


To prove (13.21) and (13.31) . we get from (j2.18p . (I2.19|) . (12.221) and (12.51) that 
/* + 00 


^ \\dx[Ti,Te,tpi:1pe]f 


dt 


c+oo 


=2 


r 


r 

p-\-oo 

d 


+ 


dt 1 


J K-i- 


J K-|_ 

Jo 


dt 


<C + C / \\dx[ipi,ipe,lfi,')pe,E,dj:[-lpi,'lpe]]\f dt < +O0. 
Jo 

On the other hand, (12.2L . (12.21) ^. (I2.2D „ and (12.51) yield 


(3.4) 


p+oo 

fj f. 

p+co 

r 

L 


dt = 2 

Jo 

/ dtdxEdxEdx 

J ]R_^ 


dt 


c+oo 


=2 


{dtipi - dt'Pe) dxEdx 


dt < +00. 


(3.5) 


Consequently, (13.4|) . (13.51) together with (12.51) gives (13.21) and (13.3|) . Then (11.261) and (11.271) follows from 
(13.2[) . (j3.3l) and Sobolev’s inequality (12.6|) . This ends the proof of Theorem ll.il □ 
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4 Appendix 

In this appendix, we will give the following inequalities stated in Lemma Id-ll repeatedlv used in the paper. 
Lemma 4.1. (i) For any function h and {k + l)j > 2, there is a positive constant C such that, 


/ \d^{u — u^,)\^\h\'^dx < ^ Sh‘^{0,t) + \\dxh{t)\\‘^ 

J IR_I_ 

(a) For any functions f h and 2(k + l)j > 3, there is a positive constant C such that, 

[ \d'f,{u - u,f)Y\hdj;f\dx < <5||5x/(t)||^ + + \\d,,h{t)\\ 

Jr+ 

(Hi) For any 6 G [0,1]^ we have 

- u,)lloo < Ce^l + 

(iv) For any 9 € [0,1], q > 10, we have 

S 


+ Iff + l) dx < + Ce%l + ln(l + k) 


and 


and 


[ \d^f\dx < ^(1 + t)-^+ Ck{l + + Ce^ (1 + t)"^+« . 

(v) For q > lOj we have 

f \g\^dx < Ce^+t (1 + + cd{l +1)-^ 

J IR..|_ 

[ \dxf\^dx < Ce^+t(1 + + C{S + e)(l + t)“^. 

Proof, {i) Using (11.141) and the following Poincare type inequalities 

\h{x,t)\ < |/i(0,t)| + x^\\dxh{t)\\, 

for (fc + l)j > 2, we have 

[ \d^{u-u^W\h\^dx 


< / \dl;{u-{k{0,t) + x\\dxh{t)\\'^) dx 

Jr+ 

<Ck(0,t) / -=—- —dx + C\\dxh{t)\\ / -=—;-— 


dx 


5hf{0,t)+ \\dxh{t)\\'^ . 

[ii) By the Young inequality and Lemma [4.II (i). for 2{k + \)j > 3, we have 

[ \d^{u-u^)\^\hdxf\dx 

]R_[. 

r § 2 {k+i)j-i 


< 6 \\dxf{t)r+c 


-h dx 


.R+ (1 + (5:c)2(fe+ib- 

<^l|5x/(t)f+ Cd2('=+i)^-3 \6hk0,t) + \\dxh{t)f 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 
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{in) From Lemma ll.3l In'], we have 


\\d:^{n'"^ -n^),d^{u''^ -m*)||oo < Cminle, (1 + t) ^}. 


Thus we have 


\\d^{rf'^-n^),dx{u^'^-u^)\\^<Ce^{l+t) 

Here we have used the fact that if 0 < C < H and 0 < C < B, then C < for any 0 < 0 < 1. 

{iv) Using (ll.23p . (11.141) . Lemma fOl fni and LemmaS^ (***), we have 

f (\f\ + \g + dy^\)dx 

J K_|_ 

<C f {dxu{u^^ — u^,) + dxU^^ {u^, — u)}dx 

JS.+ 

=C f dx[{u^^ — u^,){u — u^,)]dx + 2C ( dxU^^{u^ — u)dx 

pt /*+oo 

=2C / dxU^'^ {u^ — u)dx2C / dxU^'^ {u^ — u)dx 

Jo Jt 


X X /*+co 

<C\\dxU^^\\oo / - —dx-\-C -^ / dxU^'^dx 

Jo 1 “h Sx \ St Jt 

lU iii(i + 6 t) + lUi 

<Ce^(l + ln(l + k) + C- 


1 + 5t 

where we have used u’’^(0,t) = m* and u ^ u^: as a: —>■ +oo. 
Similarly, we can obtain that 


\dxf\dx <C f {{\dlu\ + {dxu)‘^){u'"^ - uAj + d^ud^vJ'^ + \dlu'"^\ + {d^u'"^)'^] dx 

J K_u 


<c\\dxv^^\\o. I x{\dlu\ + {dxu)^)dx + c\\dxvr-\\^\\dxu\\L^ + cuaXlUi + c\\dxvr-\\^ 


<C5(1 + t)-^ + Ce\l ++ Ce® (1 +1 ) 


-i+i 


where we have used the fact that u'^^{0,t) = u* which yields u^'^{x,t) — u^, < a;||9a;u’’^||oo. 

(v) Noticing (ll.23p and the fact that t) — u* < ||oo, and applying Lemma lOl and (I1.14p . 

we obtain that 


[ \g\^dx <C f - U*)P + - M)p}da: 

</M._j_ 

+ C\\dxU^^\\l^ f \dxu\'^x'^dx + C\\dxU^^Wl^ f {u^-u)\‘^dx 

J M,-\- JM,jf- 

<C'e^+i(l + + ^^(l + t)-^ 


and 

f Ifxl'^dx <C ( {{\dlu\^ + {dxu)^)\{u'"'" - ii*)p + \dxu\^{dxvJ"^)‘^ + \dlvd"^\^ + {Oxu"^^)^} dx 

Jr+ Jr+ 

<C\\dlvr^\\^ + C\\dxvJ'^\\lJ\dxvr^\\L^ + C'|15,^u’'^||^ f [{\dlu\‘^ + {dxu)^)x‘^ + dx 

d M+ 

^Ce^+i (1 + + C{6 + e)(l + t)-^. 


□ 
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